library(nnet)
# logit
mod1 <- glm(vs ~ hp + mpg, data = mtcars, family = binomial)
# Multinomial logit
mod2 <- multinom(factor(gear) ~ hp + am, data = mtcars, trace = FALSE)
pred1 <- predict(mod1, se.fit = TRUE, type = "response")
pred2 <- predict(mod2, se.fit = TRUE, type = "probs")